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Abstract 

Dynamical models of prototype gravastars made of anisotropic dark energy are constructed, in 
which an infinitely thin spherical shell of a perfect fluid with the equation of state p = (1 — 7)0" 
divides the whole spacetime into two regions, the internal region filled with a dark energy fluid, and 
the external Schwarzschild region. The models represent "bounded excursion" stable gravastars, 
where the thin shell is oscillating between two finite radii, while in other cases they collapse 
until the formation of black holes. Here we show, for the first time in the literature, a model of 
gravastar and formation of black hole with both interior and thin shell constituted exclusively of 



X 

dark energy. Besides, the sign of the parameter of anisotropy (pt — p r ) seems to be relevant to the 
gravastar formation. The formation is favored when the tangential pressure is greater than the 
radial pressure, at least in the neighborhood of the isotropic case (to = —1). 
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I. INTRODUCTION 



Gravastar was proposed as an alternative to black holes. The pioneer of Mazur and 
Mottola (MM) [1], consists of five layers: an internal core < r < r\ , described by the 
de Sitter universe, an intermediate thin layer of stiff fluid r\ < r < r 2 , an external region 
r > r 2 , described by the Schwarzschild solution, and two infinitely thin shells, appearing, 
respectively, on the hypersurfaces r — r\ and r = r 2 . The intermediate layer is constructed in 
such way that r\ is inner than the de Sitter horizon, while r 2 is outer than the Schwarzschild 
horizon, eliminating the apparent horizon. Configurations with a de Sitter interior have 
long history which we can find, for example, in the work of Dymnikova and Galaktionov 
After this work, Visser and Wiltshire {3] pointed out that there are two different types of 
stable gravastars which are stable gravastars and "bounded excursion" gravastars. In the 
spherically symmetric case, the motion of the surface of the gravastar can be written in the 
form (3)], 

l -R 2 + V(R) = 0, (1) 

where R denotes the radius of the star, and R = dR/dr, with r being the proper time of 
the surface. Depending on the properties of the potential V(R), the two kinds of gravastars 
are defined as follows. 

Stable gravastars: In this case, there must exist a radius ao such that 

V(R )=0, V'(R ) = 0, V"(R )>0, (2) 

where a prime denotes the ordinary differentiation with respect to the indicated argument. 
If and only if there exists such a radius Rq for which the above conditions are satisfied, the 
model is said to be stable. Among other things, VW found that there are many equations 
of state for which the gravastar configurations are stable, while others are not {3]. Carter 
studied the same problem and found new equations of state for which the gravastars are 

n n n 

stable 4| , while De Benedictis et al [5| and Chirenti and Rezzolla [6[ investigated the stability 
of the original model of Mazur and Mottola against axial-perturbations, and found that 
gravastars are stable to these perturbations too. Chirenti and Rezzolla also showed that 
their quasi-normal modes differ from those of black holes with the same mass, and thus can 
be used to discern a gravastar from a black hole. 
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"Bounded excursion" gravastars: As VW noticed, there is a less stringent notion of 
stability, the so-called "bounded excursion" models, in which there exist two radii a± and a 2 
such that 

V(R!) = 0, V'ORi)<0, V(R 2 )=0, V'(R 2 )>0, (3) 
with V(R) < for R £ (R u R 2 ), where R 2 > Ru 



Lately, we studied both types of gravastars [8|-[12|, and found that, such configurations 
can indeed be constructed, although the region for the formation of them is very small in 
comparison to that of black holes. 

Based on the discussions about the gravastar picture some authors have proposed al- 



ternative models 



13_|. Among them, we can find a Chaplygin dark star jl^j], a gravastar 



supported by non-linear electrodynamics [15] , a gravastar with continuous anisotropic pres- 
sure [16 1 and recently, Dzhunushaliev et al. worked on spherically symmetric configurations 



of a phantom scalar field and they found something like a gravastar but it was unstable 
In addition, Lobo 18] studied two models for a dark energy fluid. One of them describes a 



homogeneous energy density and the other one describes an ad-hoc monotonically decreas- 
ing energy density, although both of them are with anisotropic pressure. In order to match 
an exterior Schwarzschild spacetime he introduced a thin shell between the interior and the 
exterior spacetimes. 

Since in the study of the evolution of gravastar there is a possibility of black hole for- 
mation, we can find some works considering the hypothesis of dark energy black hole. In 
particular, Debnath and Chakraborty [7] studied the collapse of a spherical cloud, consisting 
of both dark matter and dark energy, in a form of modified Chaplygin gas. They have found 
when the collapsing fluid is only formed by dark energy, the final stage is always a black 
hole. On the other hand, Cai and Wang 22] studying the collapse of a spherically symmetric 
star, made of a perfect fluid, have found that if the fluid would correspond to a dark energy, 
black hole would never be formed. 

Here we are interested in the study of a gravastar model whose interior consists of an 
anisotropic dark energy fluid 18| which admits the isotropy as a particular case, and analyze 
the effect of the anisotropy of the pressure in the evolution of gravastars. We shall first 
construct three-layer dynamical models, and then show that both types of gravastars and 
black holes exist for various situations. We have also shown a model of gravastar and 
formation of black hole with both interior and thin shell constituted exclusively of dark 



energy. The rest of the paper is organized as follows: In Sec. II we present the metrics of 
the interior and exterior spacetimes, and write down the motion of the thin shell in the form 
of equation (TTJ). In Sec. Ill we show the definitions of dark and phantom energy, for the 
interior fluid and for the shell. From Sec. IV we discuss the formation of black holes and 
gravastars for anisotropic and isotropic fluids of different kinds of energy. Finally, in Sec. V 
we present our conclusions. 



II. DYNAMICAL THREE-LAYER PROTOTYPE GRAVASTARS 



The interior fluid is made of an anisotropic dark energy fluid with a metric given by the 



first Lobo's model 



18| 



ds 2 



-fidt 2 + f 2 dr 2 + r 2 dn 2 



(4) 



where dVL 2 = d9 2 + sin 2 (9)d(j) 2 , and 



fi 

h 



[1 - 2ar z 
1 



l+3q 



(5) 



1 - 2ar 2 ' 

where a; is a constant, and its physical meaning can be seen from the following equation 
(JTj). Since the mass is given by m(r) = 47rpo^ 3 /3 and a = ATcpo/3 then we have that a > 0, 
where po is the homogeneous energy density. Note that there is a horizon at — 1/ \/2a, 
thus the radial coordinate must obey r < r^. The corresponding energy density p, radial 
and tangential pressures p r and p t are given, respectively, by 

P = Po = constant, 

Pr = LOpo, 



Pt 



1 



4tt (1 +u)(l + 3uj)p r 2 
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CO 



(1 



(6) 
(7) 



when oj = — 1 and u = -1/3 we obtain an interior isotropic pressure fluid. 
The exterior spacetime is given by the Schwarzschild metric 

ds 2 e = -fdv 2 + f^dr 2 + rW, 

where f = 1 — 2m/ r. The metric of the hypersurface on the shell is given by 



(8) 



dsl = -dr 2 + R 2 (r)dtt 2 , 



(9) 



where r is the proper time. 

Since ds 2 = ds 2 = c?s|, we find that = = R, and 



fit - h& = 1, 

/« - T = l, 



(10) 

(11) 



where the dot denotes the ordinary differentiation with respect to the proper time. On the 
other hand, the interior and exterior normal vectors to the thin shell are given by 



n. 



-R,i,0,0), 



< = (-R,v,0,0). 



(12) 



Then, the interior and exterior extrinsic curvatures are given by 



K 

T 



K e 

TT 



K e 



-(1 - 2aR 2 )-^ +1 ^ 2 { [6(1 - 2aR 2 )^ +1 ^ 2 R 2 u + 6aR 2 i 2 u + 2aR 2 i 2 - 3t 2 u - i 2 

xaRi - (1 - 2aR 2 ) {3uj+1)/2 (-l + 2aR 2 )(Rt - Ri)} (-1 + 2aR 2 )~ 1 , (13) 

i(l-2aR 2 )R, (14) 

^sin 2 (0), (15) 

v(4m 2 v 2 - 4mRv 2 - 3R 2 R 2 + R 2 v 2 ){2m - R^mR' 3 + Rv - Rv, (16) 

-v(2m-R), (17) 

K e ee sm 2 (6). (18) 



Since [l9] 



[K ee ] = K* ee - K\ e = —M, 



where M is the mass of the shell, we find that 



M = v(2m -R) + i(l - 2aR 2 )R. 



Then, substituting equations ( TTUj) and (TTTj) into (|20|) we get 



(19) 



(20) 



M 



- Ril - 2 ~w +R2 



i/2 ( 1 - 2aR 2 + R 2 
+ R 



1/2 



(21) 



;i -2aR 2 )-^+ 2 )/ 2 ' 

In order to keep the ideas of MM as much as possible, we consider the thin shell as consisting 
of a fluid with the equation of state, p = (1 — 7)0", where a and p denote, respectively, the 
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surface energy density and pressure of the shell and 7 is a constant. Then, the equation of 



motion of the shell is given by 



IS I 



M + 8nRRp = AirR 2 [T a pu a n p ] = ixR 2 - T^ufrif) , (22) 

where u a is the four- velocity. Since the interior fluid is made of an anisotropic fluid and the 
exterior is vacuum, we get 

M + 8ttRR{1 - 7)0- = 0. (23) 
Recall that a = M/(4ttR 2 ), we find that equation (|23|) has the solution 

M = kR 2 ^~ 1 \ (24) 

where k is an integration constant. Substituting equation fpM|) into equation (|2"T|) . and 
rescaling m, b and R as, 

1 

m — y mk 2 t- 3 , 
2 

a — > ak 2 ^- 3 , 

R ->• Rk~^h, (25) 
we find that it can be written in the form of equation (CQ) 

V(R, m, co, a, 7) = ^ _ {R 2 b 2 - R 2 + 2aR* + b 2 R^~ A + b 2 R 2 ^ 2 [-2b 2 R 2 ^~ 2 - 

2V2b 2 mR - 2b 2 aR A - 2mR + 2aR 4 + b 2 R^~ 4 {b 2 - l)" 1 - 2b 2 mR^ , 

(26) 

where 

b= (1 -2ai? 2 )- (1+3w/2) . (27) 

Clearly, for any given constants m, co, a and 7, equation ( 126]) uniquely determines the collapse 
of the prototype gravastar. Depending on the initial value Rq, the collapse can form either 
a black hole, a gravastar or a spacetime filled with anisotropic dark energy fluid. In the 
last case, the thin shell first collapses to a finite non-zero minimal radius and then expands 
to infinity. To guarantee that initially the spacetime does not have any kind of horizons, 
cosmological or event, we must restrict Rq to the range, 



2m < R < 1/V2a, (28) 



where Ro is the initial collapse radius. 

Solving V(R, m, u, a, 7) = with respect to m, we can find the critical mass given by 



b 2 R 



1 



-R 2 b 2 - -b 2 R^ + aR 4 - -R 2 - V-2ab 2 R 4 f + b 2 R^~ 2 
2 2 2 



(29) 



Since the potential ( 126]) is so complicated, it is too difficult to study it analytically. 
Instead, in the following we shall study it numerically. 



III. CLASSIFICATIONS OF MATTER, DARK ENERGY, AND PHANTOM EN- 
ERGY FOR ANISOTROPIC FLUIDS 



Recently [21], an explicit classification of matter, dark and phantom energy for an 
anisotropic fluid was given in terms of the energy conditions. Such a classification is neces- 
sary for systems where anisotropy is important, and the pressure components may play very 
important roles and can have quite different contributions. In this paper, we will consider 
this classification to study the collapse of the dynamical prototype gravastars, constructed 
in the last section. In particular, we define dark energy as a fluid which violates the strong 
energy condition (SEC). From the Raychaudhuri equation, we can see that such defined dark 
energy always exerts divergent forces on time-like or null geodesies. On the other hand, we 
define phantom energy as a fluid that violates at least one of the null energy conditions 
(NEC's). We shall further distinguish phantom energy that satisfies the SEC from that 
which does not satisfy the SEC. We call the former attractive phantom energy, and the 
latter repulsive phantom energy. Such a classification is summarized in Table I. 

For the sake of completeness, in Table II we apply it to the matter field located on the 
thin shell, while in Table III we combine all the results of Tables I and II, and present all 
the possibilities found. 

In order to consider the equations (jl]) and ([7]) for describing dark energy stars we must 
analyze carefully the ranges of the parameter oj that in fact furnish the expected fluids. It 
can be shown that the condition p + p r > is violated for u> < — 1 and fulfilled for ui > —1, 
for any values of R and a. The conditions p + p t > and p + p r + 2p t > are satisfied for 
uj < — 1 and —1/3 < u> < 0, for any values of R and a. For the other intervals of oj the 



energy conditions depend on very complicated relations of R and a. See 21|. This provides 



an explicit example, in which the definition of dark energy must be dealed with great care. 
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TABLE I: This table summarizes the classification of the interior matter field, based on the energy 
conditions j^ol ]. where we assume that p > 0. 



Matter 


Condition 1 


Condition 2 


Condition 3 


Normal Matter 


p + p r + 2p t > 


P + Pr > 


p + Pt>0 


Dark Energy 


p + p r + 2p t < 


P + Pr > 


p + p t >0 






p + p r < 


p + Pt > u 


Repulsive Phantom Energy 


p + p r + 2 Pt < 


P + Pr > 


p + ft < o 






p + p r < 


p + ft < o 






p + p r < 


p + Pt > o 


Attractive Phantom Energy 


p + p r + 2p 4 > 


P + Pr > 


P + Pt<0 






p + p r < 


P + Pt<0 



TABLE II: This table summarizes the classification of matter on the thin shell, based on the 
energy conditions 20(]. The last column indicates the particular values of the parameter 7, where 



we assume that p > 0. 



Matter Condition 1 Condition 2 7 

Normal Matter cr + 2p>0 a+p>0 -lorO 

Dark Energy a + 2p < a + p>0 7/4 

Repulsive Phantom Energy a + 2p < <7 + p<0 3 



Another case was provided in a previous work 2l|. In a recent paper we have considered 
several values of u in the intervals — l<w<— 1/3 and u> > for an isotropic interior model 
and we could not found any case where the interior dark energy exists, in contrast to our 
results presented in this work. 

In order to fulfill the energy condition a + 2p > of the shell and assuming that p = 
(1 — 7)0" we must have 7 < 3/2. On the other hand, in order to satisfy the condition 
o + p > 0, we obtain 7 < 2. Hereinafter, we will use only some particular values of the 
parameter 7 which are analyzed in this work. See Table II. 

In the next sections we will discuss the different types of physical systems that we can 
find in the study of the potential V(R, m, u, a, 7). 
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IV. POSSIBLE CONFIGURATIONS 



Here we can find many types of systems, depending on the combination of the constitution 
matter of the shell and core. Among them, there are formation of black holes, sta 
gr avast ars and dispersion of the shell, as it has already shown in our previous works 



8] 



12J and listed in the table III. 



As can be seen in the figures [H El El and HI the formation of the gravastar appears as the 
unique possibility, for a given choice of the physical parameters. We can see that V(R) = 
now can have one, two or three real roots, depending on the mass of the shell. For m > m c 
we have, say, R i: where Ri+i > R4. If we choose i? > R3 (f° r m = m c we have R 2 = R3), 
then the star will not be allowed in this region because the potential is greater than the 
zero. However, if we choose R\ < Ro < R2, the collapse will bounce back and forth between 
R — Ri and R = R 2 . Such a possibility is shown in these figures. This is exactly the 
so-called "bounded excursion" model mentioned in 3], and studied in some details in 8]- 
[121 ] . Of course, in a realistic situation, the star will emit both gravitational waves and 
particles, and the potential will be self- adjusted to produce a minimum at R = R s ,„,h where 
V (R = R s tatic) = = V (R = Rstatic) whereby a gravastar is finally formed 0, U[lO] . 

Moreover, the model considered now, allows us to investigate if the anisotropy has a 
relevant role in the gravastar formation. The graphics [H El E] and H] also suggest that 
gravastar is formed for isotropic as well as anisotropic fluid cores. However, it is interesting 
to note that the sign of the parameter of anisotropy (pt — p r ) seems to be relevant to this 
formation. Those figures indicate, in a first view, that gravastar formation is favored when 
the tangential pressure is greater than the radial pressure, at least in the neighborhood of 
the isotropic case (u = —1). 

Another two completely new configurations appear from the potential studied here. They 
are represented by figures [5] and [7] and corresponding to a gravastar and a black hole, 
respectively, formed uniquely by non standard energy. The former one corresponds to a 
stable system formed by core and shell made of repulsive phantom energy. It means that a 
system constituted only by dark energy is able to achieve an equilibrium state. The second 
one reinforces this last conclusion since that it shows that even a black hole can be formed 
by a system formed exclusively by dark energy. 

Thus, solving equation (|2T1) for R(t) we can integrate R{t) and obtain R{t), which are 
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TABLE III: This table summarizes all the configurations found in our analysis. The figures are 
shown only for the new systems which results were never shown before in the literature. 



Case 


Interior Energy 


Shell Energy 


Figures 


Structures 


A 


Standard 


Standard 




Black Hole 


B 


Standard 


Dark 




Black Hole 


C 


Standard 


Repulsive Phantom 




Black Hole/Dispersion 


D 


Dark 


Standard 


mm 


Gravastar 


E 


Dark 


Dark 




Black Hole/Dispersion 


F 


Dark 


Repulsive Phantom 




Black Hole/Dispersion 


G 


Repulsive Phantom 


Standard 




Gravastar 


H 


Repulsive Phantom 


Dark 


7 


Black Hole 


I 


Repulsive Phantom 


Repulsive Phantom 


® 


Gravastar 


J 


Attractive Phantom 


Standard 




None Structure 


K 


Attractive Phantom 


Dark 




Black Hole 


L 


Attractive Phantom 


Repulsive Phantom 




None Structure 



shown in the figures [6] and [HI for the cases where the core and shell are made of dark energy. 



V. CONCLUSIONS 



In this paper, we have studied the problem of the stability of gravastars by constructing 
dynamical three-layer models of VW [3], which consists of an internal anisotropic dark energy 
fluid, a dynamical infinitely thin shell of perfect fluid with the equation of state p = (1—7)0", 
and an external Schwarzschild space. 

We also have shown that the isotropy or anisotropy of the interior fluid may affect the 
gravastar formation. The formation is favored when the tangential pressure is greater than 
the radial pressure, at least in the neighborhood of the isotropic case (co = —1). See figures 
HI El El IU where u = — 1 represent the isotropic dark energy fluid and u/-l represent the 
anisotropic dark energy fluid. 

We have shown explicitly that the final output can be a black hole, a "bounded excursion" 
stable gravastar or an anisotropic dark energy spacetime, depending on the total mass m of 
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FIG. 1: The potential V(R) and the energy conditions EC1= p + p r + 2pt, EC2= p + p r and 
EC3= p + pt, for 7 = 0, w = -1, -1.0001, -1.0003, -1.0005, a = 0.37 and m c = 0.5170643256. 
Case D: to = — 1 

the system, the parameter u, the constant a, the parameter 7 and the initial position Rq of 
the dynamical shell. All the results can be summarized in Table III. An interesting result 
that we have found is that we can have gravastar and even black hole formation with an 
interior and thin shell dark energy. 
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FIG. 2: The potential V(R) and the energy conditions EC1= p + p r + 2pt, EC2= p + p r and 
EC3= p + pt, for 7 = 0, u = -1, -0.9999, -0.9997, -0.9996, a = 0.37 and m c = 0.5170643256. 
Case D: to = — 1 




FIG. 3: The potential V(R) and the energy conditions EC1= p + p r + 2p t , EC2= p + p r and 
EC3= p + ft, for 7 = 0, w = -1, -1.0001, -1.0003, -1.0005, a = 0.45 and m c = 0.5055981510. 
Case D: to = — 1 
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FIG. 4: The potential V(R) and the energy conditions EC1= p + p r + 2p t , EC2= p + p r and 
EC3= p + jh, for 7 = -1, u = -1, -0.9999, -0.9997, -0.9996, a = 0.45 and m c = 0.5055981510. 
Case D: to = — 1 




FIG. 5: The potential V(R) and the energy conditions EC1= p + p r + 2p t , EC2= p + p r and 
EC3= p + pt, for 7 = 3, w = -3, a = 0.1 and m c = 0.0479271625. The righter vertical straight line 
denotes the radius of the horizon R^ = l/y/2a. Case I 
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FIG. 6: These figures show the dynamical evolution of collapse of a gravastar, forming a "bounded 
excursion" gravastar. We have assumed the values 7 = 3.0, to = —3.0, a = 0.1, m = 0.01 and 
initial radius Rq = R(t = 0) = 1.0. 




FIG. 7: The potential V(R) and the energy conditions EC1= p + p r + 2p t , EC2= p + p r and 
EC3= p + pt, for 7 = 7/4, oj = —3, a = 0.1 and m c = 0.1. The righter vertical straight line denotes 
the radius of the horizon Rh = l/V^a- Case H 
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FIG. 8: These figures show the dynamical evolution of collapse of a gravastar, forming a black 
hole. We have assumed the values 7 = 7/4, u = —3.0, a = 0.1, m = 0.1 and initial radius 
R = R(j = 0) = 0.9. 
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